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Scale factor duality, a truncated form of time dependent T-duality, is a symmetry of string effective action in cosmological backgrounds 
interchanging small and large scale factors. The symmetry suggests a cosmological scenario ( "pre-big-bang" ) in which two duality 
related branches, an inflationary branch and a decelerated branch are smoothly joined into one non-singular cosmology. The use of 
scale factor duality in the analysis of the higher derivative corrections to the effective action, and consequences for the nature of exit 
transition, between the inflationary and decelerated branches, are outlined. A new duality symmetry is obeyed by the lowest order 
equations for inhomogcneity perturbations which always exist on top of the homogeneous and isotropic background. In some cases 
it corresponds to a time dependent version of S-duality, interchanging weak and strong coupling and electric and magnetic degrees 
of freedom, and in most cases it corresponds to a time dependent mixture of both S-, and T-duality. The energy spectra obtained 
by using the new symmetry reproduce known results of produced particle spectra, and can provide a useful lower bound on particle 
production when our knowledge of the detailed dynamical history of the background is approximate or incomplete. 



1 Introduction 

Our starting point is the effective action (in the so-called 
"string-frame"), 
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higher orders in a'd + higher orders in 



is a "Brans-Dicke-like" scalar with uibd 



The dilaton . 
-1. 

The basic length scales of the theory are the string 
scale a' = t\ and the Planck scale e^a' = Gat = £p, re- 

lated by = g 2 strmg ~ ±a GUT (l/£ s ) = These 
relations are modified for strongly coupled string the- 
ory. We will assume that the theory is weakly coupled 
throughout the evolution. 

At early times fields may have been displaced from 
their present state, so we look for general FRW type solu- 
tions g^y = diag(— 1, a 2 (t)dxidx l ), cf> — 4>(t), and if other 
fields are present we assume that they have only time 
dependence. We obtain equations for the Hubble param- 
eter H(t) — a(t)/a(t), 4>{t), and additional fields, in par- 
ticular moduli, and matter (if present), and solve them, 
requiring that at late times the evolution has to be that 
of standard cosmology. 

In a general effective action we may represent any 
contributions in addtiinii to lowest order action by a 
"matter" LagrangianErErtil 



Seff = I ^ 
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The equations of motion arc the following 



(2) 



3H S ± X /3H 2 



s + eV Ps 



H s = ±H S ^ZH 2 + e* Ps + -e*(p s + A^C m ) (3) 
ps + 3H s (ps + Ps) = -A^Cmd, 



where 
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1 1 8C T 



2V-S 
diag(p,-p,p, -p). 



(4) 



As a result of scale factor duality (SFD) which will 
be discussed in more detail below, the solutions come in 
pairs or branches, the (+) branch vacuum (without any 
sources) satisfies 



3iJ|, 



Hs = ~7gt~{! t < to and is characterized by a future 
singularity. If the universe starts expanding according 
to the (+) branch vacuum solution, H is positive, and 
therefore its time derivative is also positive. This branch 
cannot connect smoothly to radiation dominated (RD) 
FRW with constant dilaton, <j) = 3H S + ^/3i?| + e^p s 
<t> > 0. The (— ) branch vacuum satisfies 



H S 



H 



3HI 



Hs = -j-zj^t > to an d is characterized by a past 

singularity. This branch can connect smoothly to RD 
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FRW with constant dilaton 



= 3H S - 



3fl| + e^ps 



if QHl = e^ps. 



A cosmological scenario ( "pre-big-bang" ) in which 
two duality related branches, an inflationary branch and 
a decelerated branch are smoothly joined into one cos- 
mology has been proposed Era. In this scenario the uni- 
verse quickly becomes homogeneous, isotropic, and spa- 
tially flat. The transition between the inflationary and 
decelerated branches, the so-called graceful exit tran- 
sition, is expected to occur when the universe reaches 
string scale curvature. The use of scale factor duality in 
the analysis of the higher derivative corrections to the 
effective action, and consequences for the nature of exit 
transition are outlined. 

Quantum fluctuations superimpose on top of the 
smooth classical background inhomogeneity perturba- 
tions, which are then amplified by the accelerated ex- 
pansion of the universe and materialize as particles with 
specific energy spectra later ona. 

A new duality symmetry is obeyed by theJowest or- 
der equations for inhomogeneity perturbationsH. In some 
cases it corresponds to a time dependent version of T- 
duality, interchanging small and large scale factors, in 
some cases it corresponds to a time dependent version of 
S-duality, interchanging weak and strong coupling and 
electric and magnetic degrees of freedom, and in most 
cases it corresponds to a time dependent mixture of both. 
As in other applications, duality turns out to be a power- 
ful tool for obtaining results that are inaccessible other- 
wise. In particular, lower bounds on the energy density 
of the produced particles. The energy spectrum obtained 
by truncating the solutions of the perturbation equations 
to the constant modes, is characterized by a residual du- 
ality symmetry, reproduces known results of produced 
particle spectra, and can provide a useful lower bound on 
particle production when our knowledge of the detailed 
dynamical history of the background is approximate or 
incomplete. 

2 Scale Factor Duality 

2.1 Lowest order 

To introduce SFD we look at the lowest order 4-d effective 
dilaton-gravity action 
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Integration by parts leads to 
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where ef> = <fi — 3 In a, and we have set the lapse function 
n{t) to unity. 

The action (JsJ) is invariant under the symmetry 
transformation 



«(*)- 

and H 2 are invariant, 



l/o(t), 

tfi — 61na(i), 



4>{t) -> 0(t), 

H(t) -> -H(t), 

and the equations of motion are covariant. The two 
branches describing an expanding universe are related 
to each other SFDx Time reversal. In general, for more 
complicated cosmological backgrounds the symmetry of 
the action is more complicated 13. 

2.2 Leading corrections 

The transition between the two duality related branches 
is called the graceful exit transition. It is known that 
to lowest order the two branches are separated by a sin- 
gularityLI. The emerging scenario for the exit transition 
requires classical a' corrections which can bound the cur- 
vature below the string scale, as well as quantum cor- 
rections EM. The leading classical corrections determine 
whether the solution can reach a "good" region in 0, H 
phase space. A model for the exit transition has been 
presented and therefore we know that a transition is 
possible. The question is whether string theory actually 
determines the coefficientsjsuch that a transition occurs. 

We have investigated^ effective classical corrections, 
and demanded that the action will really be an effective 
dilaton-gravity action, without additional new degrees 
of freedom EiJ. This means that we have to use actions 
that produce equations without higher derivatives. Field 
redefinitions can change that but, it is better to use a 
"frame" with no higher derivatives ensuring numerical 
stability and control. We also require an action that re- 
produces whatever string theoretic information available 
such as scattering amplitudes, perturbative beta-function 
calculations etc. As we will see this is not enough to ob- 
tain the full corrected action. If scale factor duality could 
be imposed in a practical way it would help, however, the 
situation is more complicated. 

The action including leading classical corrections is 
given by 

Slcc = ! d A x^ge-*{R + (V0) 2 + 



-[A(V^) 4 + BR 2 GB - 
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C [R^ - LgTR) V^V^ + DV 2 (^) (V^) 2 ]}. 
In covariant variables, it takes the following form 



Slcc — e 



$ f 3ff(t) a 

\ n(t) n(t) 

3 {27A + 8B + 9C + 27D) H(tf 

2 n(i) 3 

(54A + 4.B + 9C + 45L>) H(tf 4>(t) 

3 (18 A + C + 12 L>) i? (<) 2 0(t) 2 
2n(t) 3 

3 (2A + D) H(t)${t) 3 
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(3A + D) 4>(t) 4 
6n(tf 



where we have performed integration by parts to get rid 
of the h'(t), 4>"(t) and n'(t). This is possible in general 
due to the 'no higher derivatives' condition. 

Perturbative string calculation can provide two of 
the coefficients, one additional coefficient sets the over- 
all scale at which the leading corrections kick in, so one 
coefficient remains unknown. If some symmetry princi- 
ple, such as SFD could be used, the leading corrections 
action would be determined completely. For example, to 
impose naive SFD we would have to eliminate the odd 
parts of the action. We can set, for example, 
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the existence of an exact conformal field theory solution 
corresponding to the algebraic fixed point. In general 
this requires working with 2-d conformal field theories 
rather than with effective actions. I outline here som£_pf 
the possibilities to achieve progress in that directionE3. 

For highly symmetric backgrounds, such as the 
linear-dilaton deSitter background, it is possible to use 
the isometries of the background to impose additional 
symmetries on world-sheet operators, and constrain the 
beta-function coefficients. Another possibility is to use 
conformal perturbation theory to add (1,1) operators to 
an established conforrpal field theory, a linear dilaton 
flat-space background Yet another possibility is to 
start with known exact solutions in higher dimensions 
and compactify down to 4-dli2l. 



3 "S" -Duality 

The quadratic action for perturbations of any tensor field 
expanded around a cosmological dilaton-gravity back- 
ground is given by 

S pe rt = \\ d 3 xdr,S(ri) [V' 2 - (VV;) 2 ] . (7) 

The prefactor S(rj) is given by a 2m e e ^, where m and I 
depend on the type of field. For example, gravitons, dila- 
tons and moduli have m = 1, 1 = — 1, model independent 
axions have m = 1 , i ~ 1 , while Ramond-Ramond axions 
have m = 1, t = 0, and so on. We would like to compute 
the evolution of perturbation and eventually compute an 
important physical observable quantity: the spectrum of 
produced particles at late times. 



C = AA 
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However— this is only possible in a homogeneous 
background E3. In an inhomogeneous background we 
get more equations and the only consistent solution is 
A = B = C = D = Q. The conclusion is that, at the very 
least, the SFD transformation has to be modified at this 
order, making it less useful for determination of the one 
remaining coefficient. If we insist on naive SF-D in the 
homogeneous case, then it is possible to show t2l that if 
a stable algebraic fixed point exists, another non-stable 
fixed point will also exist, and that the generic solution 
will encounter the unstable fixed point first and run into 
a singularity. 

2.3 All orders 

As we have seen, additional input is required to deter- 
mine the behaviour of solutions when classical stringy 
corrections are included. The best would be to establish 



3.1 "S" -Duality symmetry 

To discuss duality symmetry of the action (0) it is more 
convenient to to use the Hamiltonian formalism. The 
Hamiltonian density corresponding to (|7j) is given by 



H 



(8) 



where the momentum conjugate to ^ is given by 

IT = S 1 *'. (9) 
The Hamilton equations of motion are first order 



n' = -^ = 5v 2 * 



(10) 
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and lead to second order equations 

n" - ^n' - v 2 n = o 
s' 

- V 2 * = 0. 

o 



(11) 

(12) 



The second equation ( |1 2| ) is commonly used in analysis 
of perturbation spectra. 

In Fourier space the Hamiltonian density is given by 

H = \J dJs-^Tl^ + Sk 2 *^^, (13) 

and the equations of motion are given by 
IT- = -STc 2 * 



9'r = S- 1 IL e . 



(14) 



"S" -duality exchanges the variables and momenta 
and at the same time sends S to its inverse, 



S ^ 5* — - S j 



(15) 



leaving the Hamiltonian, equations of motion and Pois- 
son brackets invariant. 

We are interested in a situation in which the initial 
conditions correspond to zero-point vacuum fluctuations 
of the field \&. and therefore 



(5 _1 n 2 ) = (S(W) 2 



(16) 



where (• • •) denotes ensemble average. 

The duality ( |l5| ) , contains strong- weak coupling du- 
ality as a special case. For perturbative heterotic 4-d 
gauge bosons the function S is given simply by S(rf) = 
e 0(»7) Recall that = g string , so the transformation 
S —> S = S is, at each time 77, simply the celebrated 
strong- weak coupling duality goring -* 9 string^ wnicn ap- 
pears as a part of the SL(2, Z) group, usually called S- 
duality. The transformation (Q) exchanges in this case 
electric and magnetic degrees of freedom. 



3.2 Approximate solutions 

To construct approximate solutions define n, whose 
Fourier modes are given by 



n k = s- 1 ' 2 n fc . 



(17) 



The new variables have simple transformation law under 
"S" -duality 



k^ k -> fc* fc = -n& 

n fc -» n fc = fc* fc 
s = s~ 1 . 



(18) 



The variables \I>, n satisfy the following Schrodinger- 
like equations 

+ (fc 2 - {s l ' 2 )"s- 1 ' 2 ) $ fc = 
V - [s- 1 ' 2 )" s 1 ' 2 ) n fe = 0. (19) 



Since S(r)) ~ f] a , the potentials V* = (S 1 / 3 )" S' 1 / 2 , 
Vn = (S^ 1 / 2 )" ' S 1 / 2 , if non- vanishing, are proportional 
to l/rj 2 . For k 2 > Vn, or equivalently (fcr;) 2 > 1 
(inside the horizon), we look for WKB-like approximate 
solutions 



$k(r)) = (k 2 -V*) 1/4 e 



-i J dr)'(k 2 -Vq,) 



1/2 



n fe (ry) = k (k 2 - V n ) 1/4 e 



-i J d^'^-Vh) 1 



(20) 



The advantage of looking at solutions ( p0| ) is that they 
manifestly preserve the "S" -duality symmetry of the 
equations, because the potentials V%, Vn get inter- 
changed under S — > S -1 . 

For very large A: 2 , fc 2 3> V*, Vn solutions ( |2C| ) reduce 
to correctly normalized vacuum fluctuations 



§ k ( v ) = fe~V2 e -*k>7-H<eo 
n fc (ry) = fc +1 / 2 e -<*•»+**'£, 



(21) 



where (^Oi <fi'o are random phases, originating from the 
random initial conditions. Note that because of the ran- 
dom phases, "S" -duality holds only on the average in the 
sense of eq.@. 

For k 2 < V*, Vn, or equivalently (krj) 2 < 1 (out- 

S'de the horizon), it is possible to write "exact" solutions 
It is convenient to define the functions Tcos(5 _1 , S), 



'a 



Tcos(S-\S) = l-k JdrnS-^rn) k Jd m S( m ) + ■• ■ 

+ (-i)«+ifc 2 (»-i) njs- 1 js--- js- 1 js+--- 

n—X times 
1 

Tsin{S- 1 ,S) = k J dmS^im) (22) 

+ (-i) n+1 k 2n - 3 J s- 1 l[jsjs- 1 ~.jsjs- 1 + — 1 



n — 2 times" 



4 



in terms of which the "exact" solutions take the following 
form 

$ fe (r/) = VS{ A k Tcos(S-\S)+B k Tsm(S~)S) 1 



n fe (rj) = ^j^fc TcosiS^-^-Ak Tsin(5,5- 1 )|. (23) 

Using the relations 

[Tcos(S-\S)]' = -| Tsin(5,5- 1 ) 

[TsinCS" 1 ,,?)]' = | Tcos^S- 1 ), (24) 

and similar relations for [Tcos(5, 5 -1 )] ' and 
[Tsin(S, S 1-1 )] ', it is possible to verify explicitly 
that Ilfc in eq.(p3[) are indeed solutions of eqs.(p^). 
Formally, these solutions are valid also inside the hori- 
zon, but the functions Tcos, Tsin are not well defined 
there. 

We need to match the solutions inside and outside 
the horizon and do it such that "S" -duality is respected. 
One way of doing so is to use solutions ( |20| ) inside the 
horizon, and (23) outside the horizon and match them at 
some time near horizon exit time r) ex , for which kr\ ex ~ 1. 
Taking advantage of the phenomenon of "freezing of per- 
turbations" outside the horizon we obtain the following 
result, 



1 



**(»?) 



s, 



-1/2 

cos(kri) 

1/2 



Sex 



^ ) cos(fc?7) 



1/2 

sin(kri) 

-1/2 

^ ) sin(kr]) 



(25) 



where S re — S(f] re ). The reentry time rj re is the second 
time at which krj re ~ 1. 

3.3 Energy Spectrum 

We compute an important physical observable, the 
Hamiltonian density, 



^((fcl 2 ) 



k 2 (\y k \ 2 ) 



Using the approximate solutions we obtain 



(26) 



(27) 



It is invariant under S ex — > S~ x , S re — > S^} , and overall 
rescaling of S. Note that, for a given k, (Hk) depends 
only on S ex and S re , and not on the whole evolution. 



The spectral energy distribution, dpk/dhxk 
(k 3 /a 4 )(Hk) {to = k/a) is given by 



dp(w) _ 4 \S ex {w) S re (uj) 

d In LO [ S re [U) ) S e x (w) 

~ lu a Maxj S ' 



S re 

~q 3 q 



(28) 



It has the same invariance properties as the Hamiltonian 
density. 

From cq.(|28|) we obtain model independent lower 
bound on energy density of cosmologically produced par- 
ticles. The spectrum (28) is a sum of two terms, one be- 
ing the inverse of the other. Therefore it is not possible to 
decrease the contribution of one term without increasing 
the contribution of the other. The physical origin of this 
lower bound is indeed the uncertainty principle. Recall 
that one term originates from the contribution of the 
perturbation conjugate momentum and the other from 
the contribution of the perturbation itself. The uncer- 
tainty principle says that it is not possible to decrease 
both without limits. For specific cases, the lower bound 
may be improved using some particular properties of the 
background. 

The result ( p8| ) provides an easy and a very gen- 
eral way of computing ■ The prescription is simple. 
Once the function S(r/) is known for all times, substi- 
tute for rj ex — > fc _1 , and for rj re substitute the prop- 
erly redshifted The results obtained using this sim- 
ple method reproduce knoam fPSiilts Qbtained by explicit 
complicated calculations. ^yyiaMS 
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